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ABSTRACT 

We present Chandra and XMM-Newton observations of a small sample (11 ob- 
jects) of optically-selected Seyfert 2 galaxies, for which ASCA and BcppoSAX had sug- 
gested Compton-thick obscuration of the Active Nucleus (AGN). The main goal of this 
study is to estimate the rate of transitions between "transmission-" and "reprocessing- 
dominated" states. We discover one new transition in NGC 4939, with a possible ad- 
ditional candidate in NGC 5643. This indicates a typical occurrence rate of at least 
^0.02 years -1 . These transitions could be due to large changes of the obscuring gas 
column density, or to a transient dimming of the AGN activity, the latter scenario 
being supported by detailed analysis of the best studied events. Independently of 
the ultimate mechanism, comparison of the observed spectral dynamics with Monte- 
Carlo simulations demonstrates that the obscuring gas is largely inhomogeneous, with 
multiple absorbing components possibly spread through the whole range of distances 
from the nucleus between a fraction of parsecs up to several hundreds parsecs. As a 
by-product of this study, we report the first measurement ever of the column den- 
sity covering the AGN in NGC 3393 (N H ^ 4.4 x 10 24 cm" 2 ), and the discovery 
of soft X-ray extended emission, apparently aligned along the host galaxy main axis 
in NGC 5005. The latter object hosts most likely an historically misclassified low- 
luminosity Compton-thin AGN. 
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1 INTRODUCTION 

In X-rays, obscured AGN may be classified into Compton- 
thin and Compton-thick, according to the column of absorb- 
ing matter covering the active nucleus. The threshold corre- 
sponds to a column density Nh — o"^ 1 ~ 1.5 x 10 24 cm~ 2 . 
The fact that Compton-thick Seyfert 2s are a substantial 
fraction of the whole population of Seyfert 2 galaxies, maybe 
as high as 50% (Risaliti et al. 1999), suggests that the cov- 
ering fraction of the absorbing matter is large. If a single ab- 
sorber covers a steady-state active nucleus, the classification 
of individual objects is not expected to be time-dependent. 
A review on the observational properties of Compton-thick 
Seyfert 2 galaxies has been recently published by Comastri 
(2004). 

Bona fide Compton-thick Seyfert 2 galaxies are ob- 
served in X-rays also at energies lower than the photoelectric 
cut-off. This X-ray emission is probably due to reprocessing 



of the nuclear emission by Compton-thick matter surround- 
ing the nucleus (Matt et al. 2000), and/or by hot plasma 
in the nuclear environment (Kinkhabwala et al. 2002). We 
define hereafter reprocessing-dominated Seyfert 2 galaxies 
those, whose X-ray emission in the XMM-Newton energy 
band (E < 15 keV) is dominated by reprocessing*. The 
common wisdom so far has been to identify reprocessing- 
dominated Seyferts with Compton-thick AGN. How- 
ever, very recently transitions between "Compton-thin" and 
"Compton-thick" spectral states have been serendipitously 
discovered in a few X-ray bright Seyfert 2 galaxies (Matt 

* This definition is therefore conceptually different from ( Comp- 
ton) reflection- dominated Seyfert 2s, where the emission in the 
XMM-Newton energy band is dominated by Compton-reflection 
off the far inner side of the absorber. Nevertheless, almost all 
known "reprocessing-dominated" AGN are "Compton reflection- 
dominated" . 
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et al. 2003b, and references therein). In UGC 4203, for 
instance (Guainazzi et al. 2001; Ohno et al. 2004), an 
XMM-Newton observation detected a bright (2-10 keV flux 
~ 9 x 10~ 12 erg cm~ 2 s _1 ) AGN, with a low-energy pho- 
toelectric cutoff (corresponding to Nh — 2 x 10 23 cm -2 ). 
In ASCA observations, performed about six years earlier, 
the weaker continuum and the huge K Q fluorescent iron line 
(Equivalent Width, EW ~ 1 keV) can be instead best ex- 
plained if the spectrum is dominated by the Compton echo 
of an otherwise invisible nuclear emission. Such transitions 
have been observed in both directions, and are normally ac- 
companied by substantial changes in the observed 2-10 keV 
flux. 

This discovery stimulates some fundamental questions 
on the nature of reprocessing-dominated Seyfert 2 galaxies. 
These transitions could be due in principle to a change of the 
intervening absorption. Alternatively, Seyfert 2 X-ray spec- 
tral states dominated by reprocessing may represent phases 
of low- or totally absent activity in the life of an active nu- 
cleus, as observed, for instance, in NGC 4051 (Guainazzi 
et al. 1998), NGC 2992 (Gilli et al. 2000), and NGC 6300 
(Guainazzi 2002). In these cases, the observed transitions 
require a change by at least one order of magnitude of the 
nuclear activity level. 

Transitions between "Compton-thin" and "Compton- 
thick" spectral states have been observed in 4 Seyfert 2 
galaxies so far (see Matt et al. 2003b, and references therein), 
out of about 40 objects for which multiple X-ray spec- 
troscopic measurements are available (Bassani et al. 1999; 
Risaliti et al. 2001). However, the "parent sample" is nei- 
ther homogeneous, nor complete, being substantially biased 
toward brighter (and therefore less absorbed) objects (see 
the discussion in Risaliti et al. 1999). 

We are carrying on a XMM-Newton survey of an opti- 
cally defined and complete - albeit small - sample of Seyfert 
galaxies, classified as Compton-thick according to observa- 
tions prior to the launch of Chandra and XMM-Newton. 
The primary goal of this study is to determine the rate 
of "transmission-" (i.e. Compton-thin) to "reprocessing- 
dominated" transitions t , and their typical timescale on the 
soundest possible statistical basis. This rate might be related 
to the duty-cycle of the Active Galactic Nuclei (AGN) phe- 
nomenon, at least in the local universe, if these transitions 
are due to large changes of the overall X-ray AGN energy 
output (Matt et al. 2003b). The results of this survey are 
the main subject of this paper. 



2 THE SAMPLE 

Our objects are extracted from the sample of Risaliti et 
al. (1999), which includes nearby Seyfert 2 galaxies with 
ASCA/BcppoSAX measurements of the X-ray column den- 
sity. We have restricted our analysis to those objects, whose 
[Om] luminosity is > 10~ 13 erg s _1 . As Risaliti et al. (1999) 
discuss, this choice minimizes any bias due to incomplete- 
ness. Out of the potential fourteen members of our sam- 

T Although in this paper we will refer to "transmission-" to 
"reprocessing-dominated" transitions, we search for transitions 
in both directions 



pie, priority "A" or "B" XMM-Newton observing time has 
not been allocated to five of them (IC 2560; IRAS 07145; 
NGC 5135; IC 3639; UGC 2456). We complement the XMM- 
Newton observations with two objects observed by Chandra 
(IC 2560; NGC 5135), whose data are available in the pub- 
lic archive. The galaxies discussed in this paper are listed 
in Table 1. EPIC spectra of four of the objects listed in Ta- 
ble 1 have already been individually published: NGC 1068 
(Matt et al. 2004), the Circinus Galaxy (Molendi et al. 2003), 
NGC 5643 (Guainazzi et al. 2004a), NGC 4945 (Schurch 
ct al. 2002). The Chandra observation of IC 2560 is pre- 
sented by Iwasawa et al. (2002); the Chandra observation 
of NGC 5135 is discussed by Levenson ct al. (2002, 2004). 
For the remaining five sources (NGC 1386, NGC 3393, 
NGC 2273, NGC 5005, NGC 4939) we present here for the 
first time the results of their XMM-Newton observations. 

The average distance between the Chandra/XMM- 
Newton observation and the latest ASCA/BeppoSAX one 
of the same object is ~4.9 years. 

In this paper: energies are quoted in the source reference 
frame; errors on the count rate are at the la level; uncer- 
tainties on the spectral parameters are at the 90% confidence 
level for one interesting parameter; upper limits are as well 
at the 90% confidence level; in the calculation of the lumi- 
nosities, we adopted a Hubble constant of 70 km s _1 Mpc -1 
(Bennett ct al. 2003). Preliminary results of this study are 
discussed by Guainazzi et al. (2004b). 



3 DATA REDUCTION AND ANALYSIS 

XMM-Newton data described in this paper were reduced 
with SAS v5.4.1 (Jansen et al. 2001), using the most updated 
calibration files available at the moment the data reduction 
was performed. In this paper, only data from the EPIC cam- 
eras (MOS; Turner et al. 2001; pn, Struder et al. 2001) will 
be discussed. X-ray images are generally point-like. Devi- 
ations from point-like shapes are apparent in NGC 1068 
(Matt et al. 2004), the Circinus Galaxy (Molendi et al. 
2003), NGC 4945 (Schurch et al. 2002), NGC 5005 (Sect. 5). 
Event lists from the two MOS cameras were merged before 
accumulation of any scientific products. Single to double 
(quadruple) events were used to accumulate pn (MOS) spec- 
tra. High-background particle flares were removed, by ap- 
plying standard thresholds on the single-event, E > 10 keV, 
At = 10 s light curves: 1 counts s _1 and 0.35 counts s _1 
for each pn and MOS camera, respectively. Source spectra 
were extracted from 40" circular regions around the X-ray 
nuclear source centroid, except for NGC 5643 (Guainazzi 
et al. 2004a), where a smaller region was chosen to avoid 
a serendipitous nearby bright source. Background scientific 
products were extracted from annuli around the source for 
the MOS, and circular regions in the same or nearby chips 
for the pn, at the same height in detector coordinate as 
the source location. No significant variations in any energy 
bands has been observed during the XMM-Newton observa- 
tions presented here for the first time. Spectra were binned 
in order to oversample the intrinsic instrumental energy res- 
olution by a factor >3, and to have at least 25 counts in 
each background-subtracted spectral channel. This ensures 
that the % 2 statistics can be used to evaluate the quality 



© 0000 RAS, MNRAS 000, 000-000 



Compton-thick Seyfert 2 galaxies with XMM-Newton 3 



Table 1. Sample discussed in this paper. Njj measurements in this Table were performed prior to Chandra and XMM-Newton, and are 
listed in Risaliti et al. (1999), except for the measurement in NGC 4945, which is taken from Guainazzi ct al. (2000a) 



Object 


z 




F [oui] b 


N H a 


XMM-Newton 


Exposure time 


Time span c 






(10 20 cm" 2 ) 


(cm- 2 ) 


Obs. date 


pn/MOS or ACIS (ks) 


(years) 



NGC 1068 


0.004 


3.5 


1580 


> 10 25 


29/30-Jul-2000 


61.6/66.8 


2.5 


Circinus 


0.0015 


56 


697 


(4.3±J?) x 10 24 


6/7-Aug-2001 


70.0/76.0 


3.5 


NGC 5643 


0.004 


8.3 


69 


> 10 25 


8-Feb-2002 


7.1/9.4 


4.9 


NGC 1386 


0.003 


1.4 


66 


> 10 24 


29-Dec-2002 


13.6/17.0 


6.0 


NGC 5135 


0.014 


4.6 


61 


> 10 24 


4-Sep-2001 d 


29.3 


6.6 


NGC 3393 


0.013 


6.0 


32 


> 10 25 


5-Jul-2003 


10.9/14.2 


6.5 


NGC 2273 


0.006 


7.0 


28 


> 10 25 


5-Sep-2003 


10.0/12.6 


6.5 


NGC 5005 


0.003 


1.1 


20 


> 10 24 


13-Dec-2002 


13.1/8.8 


7.0 


NGC 4939 


0.010 


3.4 


11 


> 10 25 


03-Jan-2002 


11.5/- 


5.0 


IC 2560 


0.010 


6.5 


> 4 


> 10 24 


29/30-Oct-2000 d 


9.8 


3.9 


NGC 4945 


0.002 


15.7 


> 4 


(4.4±°-§) x 10 24 


21-Jan-2001 


19.2/22.2 


1.5 



"after Risaliti ct al. (1999); derived from ASCA or BcppoSAX observation 
ft in units of lO -13 erg cm— 2 s - 1 
c minimum distance between the ASCA/BeppoSAX and the C'/i.ancira/XMM-Newton observation 

d Chandra observation 



of the spectral fitting, pn (MOS) spectra were fitted in the 
0.35-15 keV (0.5-10 keV) spectral range. 

The residuals of fits against a power-law continuum 
modified by photoelectric absorption are shown in Fig. 1 
for all the sources presented in this paper except NGC 4945 
(Schurch ct al. 2002). Notwithstanding differences, and de- 
spite the large dynamical range in observed flux, the resid- 
uals exhibit a remarkably similar pattern. Two continuum 
components can be distinguished, joining at ~2 keV (the 
only exception being the Circinus Galaxy, whose soft X- 
ray spectrum is heavily absorbed by intervening matter in 
the plane of our Galaxy). Spectra with the best statistics 
show emission-like features in the 0.5-1.5 keV energy range 
(the exceptions being in this case IC 2560, NGC 2273, and 
NGC 4939, which have the lowest signal-to-noise soft X-ray 
spectra). Above 2 keV spectra are flat, and exhibit almost 
ubiquitously intense emission line features around 6 keV (ob- 
server's frame), the only exceptions being NGC 5005, and 
NGC 4939. The latter feature is most straightforwardly ex- 
plained as iron K Q fluorescence. These lines can be better 
appreciated in Fig. 2, where we show background-subtracted 
spectra in the energy range around the K a iron line with a 
constant linear binning of about 50 eV. 

In Sect. 4 observed spectra will be compared against 
composite "two-continuum" scenarios. In these scenarios, 
the soft X-ray spectrum can be accounted for by one of the 
possible model combinations: 

• emission from an optically thin, collisionally ionized 
plasma (mekal in Xspec, Mewe et al. 1985) with free el- 
emental abundances ( "thermal scenario" hereafter) 

• a power- law with free spectral index T S oft, plus as many 
unresolved emission lines as required according to a 90% 
confidence level F-test criterion ( "scattering scenario ") 

The hard X-ray continuum will be instead accounted 
for by one of the following models: 

• a power-law with free spectral index Yhard, covered by 
photoelectric absorption with column density Nh ( "trans- 
mission scenario") 

• a "bare" {i.e. unabsorbed) Compton-reflection spec- 



trum (pexrav in Xspec; Magdziarz & Zdziarski 1995) with 
solar abundances {"(Compton-) reflection scenario") 

These simple parameterizations yield adequate fits for 
all the spectra presented in this paper. One should, how- 
ever, be aware of possible limitations inherent to this simple 
approach. High-resolution spectroscopy of nearby Seyfert 2 
galaxies (among which NGC 1068; Kinkhabwala et al. 2002; 
Brinkman et al. 2002) has convincingly demonstrated that 
soft X-ray emission is dominated - at least in some cases 
- by emission lines, with negligible contribution by an un- 
derlying continuum. Blending of these emission lines in the 
EPIC spectra can mimic a continuum emission. This point 
is discussed in larger extent by Iwasawa et al (2002). As our 
primary concern in this paper is the characterization of the 
nuclear absorber, the uncertainties induced by a purely phe- 
nomenological modeling of the soft X-ray spectrum will not 
substantially affect the core results of our paper (Guainazzi 
et al. 2004a). In the above modeling, we exclude moreover 
the possibility that the reprocessed component dominating 
the hard X-ray spectrum in the "reflection scenario" is in 
turn absorbed - e.g. by the near side outer rim or atmosphere 
of the same matter, responsible for reprocessing. This pos- 
sibility is discussed by Guainazzi et al. (2004a) with respect 
to the NGC 5643 case. In none of the other sources dis- 
cussed in this paper we have found convincing evidence for 
this possibility. However, statistics is often not good enough 
to strictly rule it out. 



4 XMM-NEWTON/ CHANDRA RESULTS 

In this Section we summarize the results of the XMM- 
Newton and Chandra (IC 2560 and NGC 5135) observations 
of the targets listed in Table 1. 

4.1 NGC 1068 

NGC 1068 is one of the X-ray brightest and best stud- 
ied Compton-thick Seyfert 2 galaxies. Its Compton-thick 
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Figure 1. Spectra (upper panels) and residuals against a power-law continuum modified by photoelectric absorption (lower panels) for 
the XMM-Newton and Chandra observations of our sample. Readers are referred to Schurch et al. (2002) for NGC 4945. Crosses: MOS: 
diamonds: pn; dots: ACIS 
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Figure 2. Background-subtracted, linearly-rebinned pn and ACIS spectra in the 5.25-7.25 keV energy range for the same objects as in 
Fig. 1 
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nature had been suggested by the prominent and multi- 
component K Q emission line complex observed by ASCA 
(Ueno et al. 1994; Iwasawa et al. 1997), and finally confirmed 
by BeppoSAX (Matt et al. 1997a). The column density of 
the absorber covering the active nucleus probably exceeds 
10 25 cm" 2 (Matt et al. 1997a). The soft X-rays are domi- 
nated by line emission following photoionization and pho- 
toexcitation by the active nucleus emission (Kinkhabwala 
et al. 2002), with little contribution from the circumnuclear 
starburst (Wilson et al. 1992). 

The EPIC spectrum of the XMM-Newton observation 
is discussed by Matt et al. (2004). Several Fe and Ni emis- 
sion lines allowed them to study in details the nature of the 
reflecting matter. Detection of iron K a Compton-shoulder 
confirms that the neutral reflector is Compton-thick. It is 
likely to be the far side inner wall of the absorber. Iron 
(nickel) overabundance of a factor about 2 (4), for lower Z 
elements when compared to solar values was measured as 
well. 



4.2 The Circinus Galaxy 

The Circinus Galaxy hosts the closest known active nucleus. 
ASCA unveiled a reprocessing-dominated spectrum (Matt et 
al. 1996). Detection of the nuclear emission in the PDS in- 
strument on-board BeppoSAX (Matt et al. 1999) allowed to 
precisely measure the column density of the absorber cover- 
ing the nucleus (Nh — 4 x 10 24 cm -2 ). In hard X-rays the 
nuclear emission is dominated by an unresolved bright core 
on scales < 8 pc (Sambruna et al. 2001). The EPIC hard X- 
ray spectra are discussed by Molendi et al. (2003) . Again, the 
measurement of iron K a Compton-shoulder - previously dis- 
covered by Chandra (Bianchi et al. 2002) - allowed them to 
identify matter responsible for the Compton-reflection dom- 
inating below 10 keV with the Compton-thick absorber 

4.3 NGC 5643 

Maiolino et al. (1998) classified NGC 5643 as a Compton- 
thick (N H > 10 25 cm" 2 ) Scyfcrt 2 galaxy, whose 0.1-10 keV 
spectrum is dominated by free electron scattering. However, 
in a later XMM-Newton pointing Guainazzi et al. (2004a) 
measured a line-of-sight column density in this object, com- 
prised between 0.6 and 1.0 x 10 24 cm" 2 . The absorber may 
be directly covering the nuclear emission or its Compton- 
reflection. Comparison with previous BeppoSAX and ASCA 
observations unveiled dramatic changes in the 1-10 keV 
spectral shape, which can be parameterized as an observed 
photon index dynamical range Ar ~ 2.0 accompanying a 
variation of the 2-10 keV flux by a factor >10. The extreme 
variability observed in the nuclear emission of this object 
indicates the revival of an AGN which was "switched-off ' 
during the BeppoSAX observation. The interpretation of 
this large variation is, however, complicated by the fact that 
the large ASCA and BeppoSAX apertures (~3') encompass 
a bright serendipitous source (christened "NGC 5643 X-l" 
by Guainazzi et al. 2004a) , apparently located in one of the 
wide spiral arm of this face-on galaxy. Understanding the 
spectral dynamics associated with the flux changes requires 
instruments capable of distinguishing the contribution of the 
two bright X-ray sources. 



4.4 NGC 5135 

We have reanalyzed the Chandra observation already dis- 
cussed by Levenson et al. (2004). Our results are substan- 
tially coincident with theirs. The ACIS-S3 spectrum is best- 
fit in the "thermal+reflection" scenario. The soft X-ray 
spectrum requires two thermal components with kT ~ 80 
and ~ 390 eV, plus an additional emission line with cen- 
troid energy E c ~ 1.78 keV. Above 2 keV the spectrum 
is Compton-reflection dominated, consistent with the AGN 
being obscured by a column density Nh i5 9 x 10 23 cm" 2 
(for an intrinsic photon index of 1.5 and a reflection frac- 
tion <0.5). The intensity of the K a fluorescent emission line 
is (5.2±2i6) x 10" 6 photons cm" 2 s" 1 , corresponding to an 
EW against the reflection continuum of 1.7±{|;| keV. The 
absorption-corrected fluxes in the 0.5-2 and 2-10 keV en- 
ergy bands are (1.9± 2 ;g) and (1.6±J;g) x 10" 13 erg cm" 2 s" 1 , 
respectively. 

4.5 NGC 1386 

The results of the XMM-Newton observation of NGC 1386 
are presented for the first time in this paper. Two of 
the baseline scenarios can be ruled out. The "scatter- 
ing+reflection" scenario can be rejected, as it produces a 
rather bad \ 2 l v = 172.7/83. The "thermal+transmission" 
scenario yields a better fit (xV" = 131.5/92). However, it 
requires a rather flat AGN spectral index (r ~ 0.5). The two 
remaining scenarios yield comparably good fits: "scatter- 
ing+transmission": x l v — 135.4/84; "thermal+reflection": 
X 2 /v = 133.4/94. In the former, the EW of the K Q iron line 
(EW ~ 1.0 keV) is too large with respect to the expected 
values for transmission through a uniform shell of material 
encompassing the continuum source (Leahy & Creighton 
1993), assuming the best-fit N H ^ 4 x 10 23 cm" 2 . In the 
latter, two thermal components are required to account for 
the bulk of the soft X-rays, alongside a Compton-reflection 
component plus iron K a iron line dominating above about 
2 keV. The best-fit parameters for the fits discussed in this 
Section are reported in Table 2. Residuals against the best- 
fit models are shown in Fig. 3. 

4.6 NGC 3393 

NGC 3393 is the object in our sample with the lowest signal- 
to-noise in the hard X-ray band. The iron line is barely 
detectable above a very weak continuum, with Ax 2 /A^ = 
6.2/1, corresponding the the 98.3% confidence level, if one 
assumes that the line is predominantly neutral.. The scat- 
tering scenario yields xt — 1-7- Thermal model for the soft 
X-ray spectra produces a significantly better fit. In the hard 
X-ray band, transmission- and reflection-dominated scenar- 
ios yield statistically comparable fits. In the former scenario 
the EW of the K Q iron line (EW = 440 ± 180) is about one 
order-of-magnitude larger than expected from the measured 
column density [N H = (7±l) x 10 22 cm" 2 , if F = 1.9]. We 
conclude therefore that Compton-reflection dominance is the 
most plausible explanation for the hard X-ray spectrum in 
this object. The lower limit on the column density covering 
the active nucleus derived from the XMM-Newton observa- 
tion [N H > 7(9) x 10 23 cm" 2 if T = 1.6 (1.9)] strictly speak- 
ing does not rule out an - albeit extreme - Compton-thin 
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Table 2. Best-fit parameters and results for the sources in Tabic 1, whose XMM-Ncwton EPIC spectral fitting results are presented for 
the first time in this paper. The legenda for the "Model" column two-letters code is as follows: the first letter indicates the scenario, 
which best accounts for the soft X-ray spectrum: scattering ("S"), or thermal emission ("T"); the second letter indicates the scenario, 
which best accounts for the hard X-ray spectrum: transmission ("T"), or reflection ("R") 







Hard X 


ray continuum 


Ed 


lission lines 






Soft X-rays 




Source 


Model 








Ia b 


EW 


fcT 


Z 










(10 23 cm- 2 ) 


(koV) 




(keV) 


(koV) 


< Z 0) 


NGC 1386 


TR 


2 5 + - 5 


> 22 


6 41± 02 
^0.03 


0.81±°;«> 


1 K + - 4 


12 + 0.05 
"■ 12± 0.02 
0.66 ± 0.03 


0.07 ± 0.02 
= Z(0.12 keV) 


133.4/94 


NGC 3393 


TR 


1.6 ± 1.2 


> 9 


6.4 C 


0.25 ± 0.14 


1.4 ± 0.8 


i 4± 0.04 


04± ' 03 
^0.02 


55.8/43 










1 84±°- 12 
^0.04 


0.15 ± 0.09 


2 2 + °- 4 


"■ ,7± 0.08 


= Z(0.14 keV) 




NGC 2273 


TR 


1.5 ± 0.4 


> is 


6.400 ± 0.010 


2 3 + 04 


0.8 ± 0.2 


< 0.06 


55.2/51 


NGC 5005 


TT 


1 6± - 7 


0.3 ±0.2 


6.4 C 


< 0.14 


< 0.24 


60 + 03 

2 °i5:9 2 

2 - J± 0.7 


30 + 14 
u.3Urt Q 30 


126.4/122 






1 5 + °- 4 








= Z(0.60 koV) 




NGC 4939 


ST 


1.5 ± 0.5 


6.4 C 


< 0.4 


< 0.07 




2.7 ± 0.4 21.1/21 








6.96 c 


< 1.1 


< 0.21 









"calculated assuming T frozen to its best-fit value for Compton-reflection dominated spectra, and a reflection fraction < 0.5 

'in units of photons cm~ 2 

c frozen 



NGC5643 NGC1386 




0.5 13 5 lO 0.5 13 5 lO 

Energy (keV) Energy (keV) 



Figure 3. Residuals in units of standard deviations against the best-fit models as in Table 2, plus NGC 5643; pn: filled dots; MOS: 
empty circles. 
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absorber. Nonetheless, its ultimate nature is confirmed by 
a reanalysis of the BeppoSAX observation (cf. Sect. 5). An 
emission line with centroid energy E c ~ 1.8 keV is required 
at the 95.1% confidence level (Ax 2 /Au = 8.4/2). This line 
may correspond to K a fluorescence of Si, which is expected 
to be produced by Compton-reflected spectra. However, its 
EW against the reflected continuum is ~5 keV, too large to 
be produced by the same Compton-reflection responsible for 
the iron emission (Matt et al. 1997b). 

4.7 NGC 2273 

For NGC 2273 the family of models where hard X-rays 
are accounted for by an absorbed power-law yield an un- 
acceptably flat intrinsic spectral index (T ~ —0.2-0.5), as 
well as an unacceptably large EW of the iron K a iron line 
(EW ~ 2.3-3.6 keV) with respect to the measured column 
density (Nh ^ 1.4-12 x 10 22 cm -2 ). Compton-reflection 
domination is a viable alternative. Modeling the soft X-rays 
with the "thermal" or the "scattering" scenario makes very 
little difference on the properties of the hard X-ray contin- 
uum or of the K a iron line, although in the latter scenario 
the photon index best-fit value is closer to standard values 
for AGN (r ~ 1.5 versus 1.2, respectively). In Table 2 we 
list the results obtained with the former. 

4.8 NGC 5005 

The XMM-Newton observation shows that the X-ray emis- 
sion is extended, and apparently elongated along a direc- 
tion close to the main axis of the host galaxy, or coinci- 
dent with an inner spiral arm, visible in the simultaneous 
OM UVW1 filter (2500-40001) image (Fig. 4). Although 
the diffuse emission is mostly associated with soft X-rays, 
the statistics is not good enough to estimate a threshold en- 
ergy, above which the X-ray emission is no longer extended. 
Assuming that the diffuse emission is due to shocked gas in 
regions of intense star formation, we have considered only 
models where at least part of the soft X-rays are due to a 
thermal component. Hard X-ray Compton-dominance is un- 
likely. A fit where the hard X-ray emission is due to a "bare" 
Compton-reflection yields a very steep intrinsic spectral in- 
dex (Yhard — 3.1). Moreover, no iron K a fluorescent line is 
detected, and the upper limit on the EW of a 6.4 keV narrow 
Gaussian profile is rather strict (< 240 eV). Transmission 
through a moderate absorber (Nh — 3 x 10 22 cm~ 2 ) is a vi- 
able alternative. The soft X-rays can be accounted for by the 
combination of two thermal components (x 2 /v = 126 .4/122) 
or of one thermal component and a scattered power-law 
(x 2 /v = 138.9/124). In Tab. 2 we show the results ob- 
tained in the former scenario. In the latter, Y^ard — 1-8, 
and Nh — 5 x 10 22 cm~ 2 . 

4.9 NGC 4939 

NGC 4939 was serendipitously located in the pn field of 
view of an observation of SAX J1305. 2-1020. Its spectrum 
is the only one in our sample, which clearly exhibits a 
soft photoelectric cut-off (cf. Fig. 1). Indeed, the "trans- 
mission" scenario accounts well for the hard X-rays, with 
N H ^ 1.5 x 10 23 cm~ 2 . Modeling the soft X-rays with a 



single thermal component (kT ~ 0.7 keV) or a steep power- 
law (F ~ 2.7) yields fits of equivalent statistical quality: 
X 2 /^ = 18.9/20, and 21.1/21, respectively. The best-fit pa- 
rameters and results for the latter are shown in Table 2. An 
emission line is additionally required at the 94.0% confidence 
level only (A^ 2 /A^ = 5.4/2). Its centroid energy is incon- 
sistent with emission from neutral iron: E c — 6.71±o;2ci keV. 
If, following Maiolino et al. (1998; cf. Sect. 5 as well), we in- 
terpret this shift of the centroid energy as due to a blend of 
neutral and a H-like transitions, the 90% upper limits on the 
EW of either component are 70 eV and 210 eV, respectively. 

4.10 NGC 4945 

The XMM-Newton observation of NGC 4945 is discussed by 
Schurch et al. (2002). The galaxy core has a complex mor- 
phology. It is dominated by reprocessing, as the nucleus is 
covered by a thick absorber (Nh — 4 x 10 24 cm -2 ; Done 
et al. 1996). Compton-reflection from the inner side of an 
edge-on torus leaves its imprinting in the hard X-ray spec- 
trum through a 1.6 keV EW K a emission line, consistent 
with previous findings (Guainazzi et al. 2000a; Madejski et 
al. 2000). In soft X-rays, multi-temperature emission from a 
nuclear starburst dominates, a two temperature model yield- 
ing kT ~0.9 keV and kT ~ 6.9 keV. The hard X-ray emission 
exhibits a resolved morphology, suggesting that part of the 
gas in the starburst region is exposed to the AGN radiation 
as well. 

4.11 IC 2560 

IC 2560 has not been observed by XMM-Newton. The 
results of a Chandra observation of this target are dis- 
cussed by Iwasawa et al. (2002). A model constituted by 
a two-component thermal emission plus a Compton reflec- 
tion dominated spectrum (with F = 1.9) is an adequate 
description of the spectrum (x 2 /v = 61.2/74). The EW 
of the K a iron line is ~ 3.5 keV. In principle, a statis- 
tically equivalent fit is obtained if the "bare" Compton- 
reflection component is substituted by an absorbed power- 
law (x 2 l v = 59.9/73). However, in this scenario the K a 
iron line EW (~ 8.0 keV) is almost two orders of mag- 
nitude too large than expected in transmission from the 
measured column density (Nh — 5 x 10 22 cm -2 ). The 
Chandra observation therefore supports the interpretation 
of the IC 2560 ASCA spectrum as hard X-ray reprocessing- 
dominated (Risaliti et al. 2000), against the interpretation of 
the same data in terms of a Compton-thin absorber covering 
the nuclear emission given by Ishihara et al. (2001) 

4.12 Fluxes and luminosities 

In Table 3 we present the observed fluxes in the 0.5-2 keV 
and 2-10 keV energy ranges for all the sources in Table 2. 
The corresponding luminosity corrected for Galactic ab- 
sorption in the soft X-ray band ranges between 10 41 to 
10 43 erg s _1 (cf. Fig. 10). In the hard X-ray band, the deter- 
mination of the intrinsic AGN luminosity is impossible for 
all cases where only lower limits on the nuclear absorbing 
column density exist. For Compton-thick sources they are 
anyhow plagued by large uncertainties. For the two sources 
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Figure 4. Left: pn intensity contours in the 0.2-15 keV energy band superposed to the OM exposure with the UVW1 filter for NGC 5005. 

// 

The pn image is smoothed with a 6 kernel boxcar function. Contours represent 15 logarithmically equispaced count levels from 0.17 to 
18.29 on the smoothed image. Right: zoom of galaxy surface in the UVW1 filter. The position of the nucleus, of the inner and the outer 
arms are labeled 



Table 3. Observed fluxes for the sources listed in Table 2. Units 



Source 



0.5-2 keV 



2-10 keV 



NGC 1386 
NGC 3393 
NGC 2273 
NGC 5005 
NGC 4939 



1.8± 



_0.9 
-0.5 
0.43±°;?| 

0.47 ±0.03 
0.12 ±0.04 



0.27 ±0.05 

0.09±°$ 
0-69±°il 
0.51 ± 0.06 

°- o;c 0.30 



which are Compton-thin in XMM-Newton observations, the 
2-10 keV luminosity is 1.8 x 10 42 erg s _1 (NGC 4939), and 
1.2 x 10 40 erg s" 1 (NGC 5005), respectively. 



Table 4) , where a delay in the response of a variable primary 
continuum probably occurs. 

In the following, some additional details are given on the 
analysis of the ASCA/BeppoSAX observations, whenever 
our analysis reaches further or different conclusions with re- 
spect to what published in the literature, or shown by the 
XMM-Newton observations. 

NGC 5135: the absorption-corrected 0.5-2 keV flux dur- 
ing the January 1995 ASCA observation was (1.2±jj;f) x 
10 -12 erg cm -2 s _1 . At face value this is one order of magni- 
tude larger than measured by Chandra 6.6 years later. How- 
ever, the difference is at the 1-a level only, if the statistical 
uncertainties are taken into account. The other spectral pa- 
rameters are consistent with the Chandra results, with large 
errors. 



5 COMPARISON WITH ASCA/BEPPOSAX 
RESULTS 

In this Section we compare the results of the Chandra and 
XMM-Newton observations described in Sect. 4 with prior 
ASCA and BeppoSAX measurements. All the spectra de- 
scribed in this Section were extracted from calibrated and 
linearized event lists available in the public archive, and re- 
analyzed by us. Whenever more than one observation was 
available for a given source, we have considered the latest (in 
no case significant spectral variability was observed across 
different ASCA/BeppoSAX observations, with the only ex- 
ception of NGC 1068 (Guainazzi et al. 2000b; Colbert et al. 
2002). This exception does not substantially affect any of 
the results discussed in this paper. 

Variability in the soft (0.5-2 keV) and hard (2-10 keV) 
X-ray flux is generally restricted to a factor of ±3 (Guainazzi 
et al. 2004b). The intensities (but not the EWs: see Sect. 6 
below) of the K a iron lines are consistent with a factor ±2 
as well. The only exception is NGC 3393 (cf. Table 2, and 



NGC 3393: NGC 3393 is one of the few targets in our sam- 
ple, which was detected by the PDS instrument on board 
BeppoSAX above 15 keV (count rate: 0.39 ± 0.09 s" 1 ). 
The BeppoSAX observation is discussed by Maiolino et al. 
(1998). They interpret the BeppoSAX spectrum as due to a 
Compton-thick source, with a column density > 10 25 cm -2 . 
The overall poor statistics of the BeppoSAX observation 
prevented them from applying more complex models. How- 
ever, the PDS data points in their Fig. 1 lay systematically 
above the extrapolation of the best-fit model in the 2-10 keV 
band. 

We have first applied the best model of the 
XMM-Newton observation to the BeppoSAX spectra. As 
NGC 3393 is undetected by the LECS instrument below 
2 keV, we kept the parameters of the thermal components 
and of the E c ~ 1.8 keV emission line frozen to the XMM- 
Newton best-fit values, as the contribution of these com- 
ponents is negligible in the MECS-PDS energy bandpass. 
The best-fit intrinsic spectral index face value turns out to 
be P ~ 0.7. Although the quality of the fit is acceptable 
(x 2 l v = 21.4/17), this flat index - still consistent with the 
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Figure 5. Spectra (upper panel) and residuals in units of stan- 
dard deviations (lower panel) when the XMM-Newton best-fit 
spectrum is applied (with T = 1.6) to the BcppoSAX spectra of 
NGC 3393 (MECS: crosses; PDS: filled dots) 



Table 4. Best-fit parameters and results when the XMM-Newton 
model (expanded with a power-law covered by an absorber of 
column density Nh to account for the emission in the PDS energy 
band) is applied to the BeppoSAX spectra of NGC 3393. Details 
in text. 



r 

N H (10 24 cm" 
Fe K a line: 
E c (keV) 
Ic a 

EW (keV) 



2.8±J; 2 

2.5 

1 



'■) 4.4± 2 



6.58±[Ji? 
1.4 ±0.8 

4±2 
13.6/15 



a in units of 10 photons cm 



XMM-Newton results within the large statistical uncertain- 
ties - is suggestive of additional spectral complexity. If F 
is fixed to the XMM-Newton best-fit value (1.6), the PDS 
counts are largely underpredicted (Fig. 5). The difference 
is even larger if more typical values > 1.9 are used. We 
conclude that the flux in the PDS band is dominated by 
the nuclear emission piercing through a Compton-thick ab- 
sorber, with Nh < 10 25 cm~ 2 . Adding an absorbed power- 
law to the XMM-Newton best-fit model yields an improve- 
ment in the quality of the fit at the 96.7% confidence level 
(Ax 2 /A^ = 7.8/2), with N H ~ 4x 10 24 cm~ 2 , and a slightly 
steeper intrinsic spectral index. The best-fit parameters are 
shown in Table 4. These results confirm that the column 
density covering the NGC 3393 nucleus is indeed Compton- 
thick, although not large enough to fully suppress the nu- 
clear emission. 

NGC 4939: NGC 4939 was classified as a Compton- 
reflection dominated Compton-thick AGN by Maiolino et al. 
(1998), on the basis of the very flat spectral index obtained 
in the "transmission-scenario" , and the fact that this model 
underpredicts the emission in the PDS band (13-200 keV 
count rate: 0.20 ± 0.06 s _1 ). We have reanalyzed the same 
data, obtaining results which are basically consistent with 
theirs. The EW of a single Gaussian profile accounting for 
the observed iron emission line (~ 750 eV) is indeed too large 



Table 5. Best-fit parameters and results for the best-fit to the 
BeppoSAX spectra of NGC 4939. Details in text. 



1 hard 


1.90±g}6 


F soft 


3.5±°| 


N H (10 24 cm- 2 )" 


> 2 


0.5-2 keV flux 6 


0-43±g;i| 


2-10 keV flux 6 


1.6 ±0.2 


Fe K a lines: 




-fc,6.4 C 


1.1 ±0.6 


EW C ,6A (eV) 


490 ± 270 


-fc,6.96 C 


1 2+ ' 6 


EW c , 6 .m (cV) 


460± 2 ™ 


x 2 ' l v 


45.4/55 



"assuming Y hard = 1.9 

6 in units of 10 — 12 erg cm -2 s _1 

c in units of 10 -5 photons cm~ 2 
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Figure 6. 3-10 keV total spectral energy distributions of 
NGC 4939 based on the 1997 BeppoSAX and 2002 XMM-Newton 
observations (solid lines). The dot-dashed lines indicate the up- 
per limits for the 6.4 keV and 6.96 keV K a fluorescent emission 
lines in the XMM-Newton observation. For display purpose only, 
we have attributed to all emission lines a formal width equal to 
the intrinsic instrumental energy resolution: 150 eV and 500 eV 
for XMM-Newton (pn) and BeppoSAX (MECS), respectively 



with respect to the measured column density in the trans- 
mission scenario (Nh — 1.3 x 10 23 cm -2 ). A fit where hard 
X-rays are dominated by a "bare" Compton-reflection is ex- 
cellent (x 2 /v — 45.4/55; Table 5). The reflection-dominated 
state is confirmed by the large EW of the neutral component 
of the iron K a fluorescent line (EW ~ 500 eV). The H-like 
component exhibits a comparable EW. The combination of 
hard X-ray continuum and iron emission line EW points to 
a transition between a "reprocessing-" and a "transmission- 
dominated" state occurring between the January 1997 Bep- 
poSAX and the March 2001 XMM-Newton observation. A 
comparison between the 3-10 keV spectral energy distribu- 
tion, based on the best-fit BeppoSAX and XMM-Newton 
models, is shown in Fig. 6. It is interesting to observe that 
the soft (E < 2 keV) X-ray flux decreased by a factor ~3.5 
between the BeppoSAX and the XMM-Newton observation. 
This supports our interpretation of the soft X-ray emission 
in this object as due to scattering of the primary nuclear 
continuum, which was mirroring a previous phase of strong 
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Figure 7. Absorbing column density Nh versus 2-10 kcV to 
0[m] flux ratio for the sample of Seyfert 2 galaxies after Bassani et 
al. (1999). The empty circle indicates the location corresponding 
to NGC 5005 for the column density observed by XMM-Newton, 
and the 0[in] fluxes reported in the literature (details in text) 



6 DISCUSSION 

6.1 How much do we know of Compton-thick 
Seyfert 2 galaxies? 

The safest way to identify a Compton-thick Seyfert 2 galaxy, 
and describe - even at the simplest phenomenological level - 
the basic X-ray properties of its nuclear emission is to detect 
the primary continuum piercing through the Compton-thick 
absorber. This requires measurement above 10 keV, which 
have been possible so far only on the ~10 objects detected 
by the PDS instrument on-board BeppoSAX. For all the 
remaining known ~40 Compton-thick Seyfert 2 (Comastri 
2004) the classification relies on indirect evidence, such as 
the flatness of the hard X-ray continuum, the EW of the K a 
iron fluorescent emission line(s), or anomalous low values of 
the ratio between the flux in the 2-10 keV energy band and 
in other wavelengths. 

Waiting for an X-ray detector of better > 10 keV sen- 
sitivity than the PDS, the robustness of the criteria used 
to identify Compton-thick objects can be tested with the 
improved sensitivity that the XMM-Newton optics offer. In 
our sample, the Compton-thick nature is confirmed for all 
objects, except NGC 5005 (N H ~ 1.5 x 10 23 cm -2 ), and - 
marginally - NGC 5643 (Guainazzi et al. 2004a; N H = 6- 
10 x 10 23 cm -2 ), apart from NGC 4939, obviously. This is a 
potentially important result, as it underlines the perspective 
to extend the search for Compton-thick objects at higher 
redshift (Fabian et al. 2002). Although classification of an 
individual object may be subject to uncertainties even when 
large EW iron lines are detected, the method is robust. 



AGN activity during the BeppoSAX observation. 

NGC 5005: NGC 5005 was declared Compton-thick by 
Risaliti et al. (1999) on the basis of the low X-ray versus 
0[lll] luminosity ratio, although no evidence for either a flat 
hard X-ray spectrum or for a K a iron line was observed 
in the ASCA spectrum. The upper limit of the EW of the 
latter feature (900 eV) was still consistent with heavy obscu- 
ration. None of the criteria adopted to classify this source 
as a Compton-thick object resists scrutiny after the XMM- 
Newton observation. The upper limit on the K a iron line 
EW is strict (< 240 eV). The application of the best-fit 
XMM-Newton model (cf. Tab. 2) to the ASCA yields a good 
fit (x 2 /f = 216.1/227), showing that the ASCA data have 
not enough statistics to distinguish a "transmission- " from 
a "reprocessing-dominated" scenario on the basis of the X- 
ray continuum shape. Literature measurements of the 0[lll] 
flux - once corrected for optical reddening using the prescrip- 
tion in Bassani et al. (1999) - span a rather large interval, 
between 0.3 and 20 x 10~ 13 erg cm~ 2 s _1 (Shuder & Os- 
terbrock 1981; Dahari & De Robertis 1988; Ho et al. 1997; 
Risaliti et al. 1999). This interval is consistent with 2-10 keV 
versus 0[lll] ratio values observed in Compton-thin as well 
as Compton-thick Seyfert 2s (Fig. 7) . We therefore conclude 
that NGC 5005 is most likely a mis-classified Compton-thin 
Seyfert 2. 



6.2 The moderately unstable temper of heavily 
obscured AGN 

The scope of this paper is comparing the X-ray spectral 
properties of a complete, unbiased sample of Compton-thick 
Seyfert 2 galaxies observed with Chandra and XMM-Newton 
with prior measurements. The main scientific goal is to es- 
timate the rate of transitions between "transmission- " and 
"reprocessing-dominated" spectral states. 

These transitions were serendipitously discovered on a 
few nearby active nuclei, once a larger database of X-ray ob- 
servations allowed us some knowledge of the X-ray history 
of a wider sample of AGN. These transitions affect Seyfert 2 
galaxies (NGC 2992; Gilli et al. 2000; NGC 1365; Iyomoto 
et al. 1997, Risaliti et al. 2000; UGC 4203, Guainazzi et al. 
2001; NGC 6300 Guainazzi 2002), as well as other AGN 
(NGC 4051; Guainazzi et al. 1998, Uttley et al. 1999; 
PG 2112+059; Gallagher et al. 2004). It has been claimed 
that variability of the absorbing column density by a factor 
50 ± 30% on timescales < 1 year is common in obscured 
AGN (Risaliti et al. 2002). In one of the best-monitored 
cases (NGC 3227; Lamer et al. 2003), the symmetric pro- 
file of the absorption light curve clearly suggests an inter- 
pretation in terms of line-of-sight crossing by an individ- 
ual cloud. The transitions we are discussing in this paper, 
however, represent a different phenomenology, whereby the 
apparent variation of the absorbing column density is of at 
least one order-of-magnitude, and the state corresponding 
to the lower X-ray flux is fully reprocessing-dominated. 

The main conclusion of this study is summarized in 



© 0000 RAS, MNRAS 000, 000-000 



12 M.Guainazzi et al. 




XMM-N/Chandra line EW (keV) 

Figure 8. EW of the K a neutral iron fluorescent emission line 
measured by ASCA/BeppoSAX as a function of the same quan- 
tity measured by CTtartdra/XMM-Newton. The locus correspond- 
ing to "static" Compton-thick objects is indicated by the shaded 
area. Filled circles: objects of our sample; open square: UGC 4203, 
the prototypical "Phoenix Galaxy" (Guainazzi et al. 2002). The 
solid (dashed) line represents the identity (ifactor 3) locus 



Fig. 8, where we use the K a iron line EW - measured by 
BeppoSAX/ASCA and C7iandra/XMM-Newton - as a hall- 
mark for bona fide Compton-thick AGN. Only "narrow" 
or "unresolved" components corresponding to neutral or 
mildly ionized iron transitions have been used in the cal- 
culation of the EWs shown in Fig. 8. We do not consider 
the contribution of Compton-shoulder, which is likely to 
be £ 20% in most cases (Matt 2002). The locus corre- 
sponding to "Compton-thick" objects is conservatively bor- 
dered by the line EW = 300 eV, corresponding to the 
brightest emission line produced by transmission through 
a Compton-thin screen, covering a 2-7T solid angle (Leahy & 
Creighton 1993). Out of the 10 Compton-thick Seyferts dis- 
cussed in this paper, we find evidence for only one transition: 
NGC 4939. This implies that a typical timescale for these 
transitions should be ~50 years, on the basis of the average 
separation between the ASCA/BeppoSAX and the closest 
CVianiira/XMM-Newton observations. Of course, given the 
size of our sample, this has to be regarded as no more than 
an order-of-magnitude estimate. However, it is consistent 
with previous determinations, based on inhomogeneous and 
incomplete samples. 

The origin of the spectral changes occurring when an 
AGN transforms its appearance from a "transmission-" to a 
"reprocessing-dominated" state is still not fully elucidated. 
In principle, high-quality, high-resolution measurements, fol- 
lowing the onset of the variability or the AGN recovery after 
a prolonged "off-state" should be decisive. 

In NGC 6300 and NGC 2992 we have the strongest ev- 
idence that transitions from transmission- to reprocessing- 
dominated states are due to a change of the optical path 



through which the nucleus is being observed, due to a tem- 
porary interruption of the nuclear activity. In NGC 6300 this 
is suggested by the very large Compton-reflection continuum 
observed by BeppoSAX (Guainazzi 2002). In NGC 2992, 
the history of the X-ray emission is comparatively well sam- 
pled. The X-ray flux decreased by about a factor of 20 from 
the first HEAO-1 detection in the early 80s, up to the 1994 
ASCA observation, just to experience a factor ~15 recovery 
by 1999 (Gilli et al. 2000). To these two cases, we might 
add the large dynamical range of the NGC 5643 AGN X- 
ray output (Guainazzi et al. 2004b), a promising transition 
candidate. Unfortunately our knowledge of the X-ray his- 
tory of NGC 4939, and UGC 4203 is too poor, for us to be 
able to draw any final conclusions on this issue. Recently 
Ohno et al. (2004) proposed a change by a factor > 5 of the 
absorber column density as the most likely explanation in 
the latter. On the other hand, multiple X-ray observations 
of optically defined samples of unobscured AGN have not 
detected a significant rate of large historical X-ray flux vari- 
ations (compare, for instance, Laor et al. 1997, George et al. 
1998, and Porquet et al. 2004). When one is dealing with X- 
ray unobscured AGN a bias toward brightest, less obscured 
object, may prevent us from discovering the fraction of X-ray 
unobscured counterparts to our "transient" Compton-thick 
Seyferts. However, this may as well suggest an alternative 
interpretation in terms of changes in the properties of line- 
of-sight gas. This scenario will be investigated in the next 
Section. 

6.3 The recovery of fossil AGN as probe of the 
circumnuclear medium 

"Transmission-" to "reprocessing-dominated" transitions 
can be used to probe some properties of the gas in the nu- 
clear environment in highly obscured AGN. The method - 
based on Monte-Carlo simulations - is described by Matt et 
al. (2003b). The application of this method to the five known 
transitions is summarized in Fig. 9, where we show the 2- 
4 keV versus 4-10 keV flux softness ratio for the Compton- 
reflection component in the reprocessing-dominated state 
against the measured column density in the transmission- 
dominated state. The solid line represents the expected cor- 
relation according to simulations. In the ASCA observation 
of NGC 2992, and possibly in UGC 4203, the spectrum ob- 
served during a reprocessing-dominated state is too hard to 
be due to Compton-reflection by matter with the same col- 
umn density as measured by the soft photoelectric cut-off 
during transmission-dominated states. A similar argument 
can be applied to NGC 6300, on the basis of its X-ray spec- 
trum above 10 keV (Guainazzi 2002; Matt et al. 2003b). The 
gas responsible for line-of-sight absorption should be there- 
fore different in density from the gas responsible for reflec- 
tion, the latter being most likely located at the far inner side 
of the molecular "torus" . This may be explained by a largely 
inhomogeneous compact (i.e., 1 pc) nuclear absorber (Ohno 
et al. 2004), or by Compton-thin absorption occurring on 
much larger scales than the nuclear "torus", e.g. in dusty 
regions associated with starburst formation (Weaver 2001), 
or with the host galaxy (Maiolino & Rieke 1995; Malkan 
et al. 1998). An extension of the Seyfert Unified scenario, 
which encompasses the latter possibility is discussed by Matt 
(2000). It is also noteworthy that in 3 out of 5 cases the SR in 
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Figure 9. 2-4 keV to 4-10 keV flux softness ratio (SR) in 
the reprocessing-dominatcd state against the Nh measured in 
transmission-dominated states for the 5 objects, where a transi- 
tion between the two states has been reported in the literature (cf. 
Sectc. 6.2), plus NGC 4939 (this paper). The solid line represents 
the expected correlation on the basis of Monte-Carlo simulations 
of Compton-refiection spectra (Matt et al. 2003b) 



Fig. 9 is softer then expected by a pure Compton-refiection 
spectrum even in the reprocessing-dominated states. This 
suggest a still not-negligible contamination by a softer com- 
ponent, e.g. a not fully "switched-off" AGN continuum. This 
possibility is confirmed by a detailed spectral analysis of the 
BeppoSAX observation of NGC 2992, the X-ray brightest 
among the objects displayed in Fig. 9, where the recovery 
of the AGN is witnessed by a comparatively small value of 
the normalization ratio between the reflected and the trans- 
mitted component [R ~ 4(fi/27r)] with respect to a bare 
Compton-refiection dominated state. 

In summary, the spectral properties of "transmission-" 
to "reprocessing-dominated" transitions indicate that ob- 
scuring matter in AGN is far from being homogeneous in 
space or time. Yet, one cannot discard the idea of a com- 
pact but inhomogeneous pc-scale "torus" (Antonucci 1993), 
or disk outflow (Elvis 2000). Such inhomogeneities might 
be the ultimate responsible for these transitions. However, 
the fact that photoionized- or starburst-dominated soft X- 
ray emission in several Seyfert 2 galaxies is absorbed as well 
(Iwasawa & Comastri 1998; Matt et al. 2001; 2003a), along- 
side our knowledge of the X-ray history in NGC 6300 and, 
above all, NGC 2992, suggests that an important contribu- 
tion to X-ray obscuration comes from matter associated to 
the host galaxy, beyond the innermost pc around the central 
engine. 



6.4 The origin of the soft X-ray emission 

Our program was not specifically tuned to investigate the 
origin of the soft X-ray emission in our sample. The an- 
swers that we get from the data on this issue are therefore 
inevitably ambiguous in almost all cases. The only excep- 
tion is NGC 5005 (the only non Compton-thick source in 
the sample). In this case, the soft X-ray emission is clearly 
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Figure 10. FIR luminosity (as defined in David et al. 1992) 
against the 0.5—4.5 keV luminosity for the Compton-thick ob- 
jects of our sample (filled dots), and the remaining Compton-thick 
Seyfert 2 galaxies of Bassani ct al. (1999). The solid line indicates 
the empirical correlation found by David et al. (1992) for starburst 
galaxies; the dashed lines indicates the 90% confidence level on 
this correlation. 'CG": the Circinus Galaxy 



extended on scales comparable with the optical size of the 
galaxy, and roughly aligned with its main axis, or with an 
inner arm UV structure. For all the other cases, the two 
proposed scenarios are equally viable on the basis of the 
EPIC data alone. Whenever high-resolution spectroscopic 
data are available, soft X-rays appear to be dominated by 
emission lines following photoionization and photoexcitation 
by the primary AGN emission (Kinkhabwala et al. 2002; 
Sambruna et al. 2001; Bianchi et al. 2001), with little or 
no contribution by nuclear starbursts. On the other hand, 
the 0.5-4.5 keV X-ray luminosities are generally consistent 
with the correlation with the Far InfraRcd (FIR) luminos- 
ity empirically determined by David et al. (1992) on a large 
sample of starburst-dominated galaxies (Turner et al. 1997; 
Maiolino et al. 1998; Fig. 10), the only discrepant objects 
being NGC 4939, and NGC 4945. This, however, holds as 
well to objects (such as the Circinus Galaxy), for which it 
is unlikely that soft X-rays are dominated by starburst. On 
the other hand, high-resolution imaging of NGC 4945 with 
Chandra shows that soft X-rays are likely to be dominated 
by thermal emission from starbursts, alongside a starburst 
mass-loaded superwind (Schurch et al. 2002). While we refer 
to Guainazzi et al. (2004a) for a more detailed discussion on 
this point, we conclude for the time being that a through- 
out determination of the physical properties of the plasma 
dominating the soft X-rays in obscured AGN requires deep 
exposures with high-resolution detector, which are currently 
possible only on a limited number of objects. 
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7 CONCLUSIONS 

The main result of this paper is the estimate of the occur- 
rence rate of "transmission-" to "reprocessing-dominated" 
state transitions, on the most unbiased and complete exist- 
ing sample of Compton-thick Seyfert 2 galaxies (cf. Fig. 8). 
We have discovered 1 (new) transition out of a sample of 10 
bona-fide Compton-thick objects, with an average time span 
between pre- and post XMM-Newton and Chandra launch 
observations of about 5 years. The statistics is still too 
small to determine anything more accurate than an order- 
of-magnitude estimate for the occurrence rate. Bearing this 
caveat in mind, it seems that once or twice every century 
we might be forced to change our X-ray absorption classifi- 
cation for each highly obscured AGN. 

With respect to the mechanism responsible for these 
transitions, our results are consistent with the discussion in 
Matt et al. (2003b). Although an explanation in terms of 
varying line-of-sight column density cannot be ruled out, 
the fact that in the best studied cases these transitions are 
associated with large (>10) fluctuations of the AGN X-ray 
output suggests that they are due to a change of the op- 
tical path through which we observe the nucleus. For the 
transition specifically discussed in this paper - NGC 4939 - 
the factor of 2 variation of the observed 2-10 keV band flux 
hides a larger variation of the AGN intrinsic power, as its 
true luminosity in the reprocessing-dominated state is un- 
known. If this is the correct interpretation, the transition 
occurrence rate translates immediately into a duty-cycle of 
the AGN phenomenon in the local universe. 

Independently of the ultimate mechanism responsible 
for these transitions, comparison of their spectral properties 
with Monte-Carlo simulations demonstrates that obscuring 
gas in absorbed AGN cannot be distributed in a space- or 
time- homogeneous structure. Again, a compact but inho- 
mogeneous "torus" cannot be ruled out. However, there is 
mounting evidence that gas in regions of intense star for- 
mation and dust in the host galaxy play a major role, and 
might be ultimately responsible for the bulk of Compton- 
thin X-ray absorption in AGN. 
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